arXiv:l504.00049vl [cond-mat.mes-hall] 31 Mar 2015 


Start 


Splitting of the monolayer out-of-plane A'^ Raman mode in few-layer WS2 


Matthias Staiger/ ’[^ Roland Gillen,^ Nils Scheuschner,^ Oliver Ochedowski,^ Felix 
Kampmann,^ Marika Schleberger,^ Christian Thomsen/ and Janina Maultzsch^ 

^Institut fur Festkorperphysik, Technische Universitdt Berlin, 

Eugene-Wigner-Bldg. EW 5-4, Hardenbergstr. 36, 10623 Berlin, Germany 
^Eakultat fiir Physik and Cenide, Universitdt Duisburg-Essen, Lotharstr. 1, 47057 Duisburg, Germany 

We present Raman measurements of mono- and few-layer WS 2 . We study the monolayer A'^ mode 
around 420 cm“^ and its evolution with the number of layers. We show that with increasing layer 
number there is an increasing number of possible vibrational patterns for the out-of-plane Raman 
mode: in N-layer WS 2 there are N F-point phonons evolving from the A'^ monolayer mode. For 
an excitation energy close to resonance with the A excitonic transition energy, we were able to 
observe all of these N components, irrespective of their Raman activity. Density functional theory 
calculations support the experimental findings and make it possible to attribute the modes to their 
respective symmetries. The findings described here are of general importance for all other phonon 
modes in WS 2 and other layered transition metal dichalcogenide systems in the few layer regime. 


PACS numbers: 78.30.-j, 78.67.-n, 63.22.Np, 71.35.Gg 


Introduction 

The last few years have seen a spectacular increase 
in interest in layered transition metal dichalcogenides 
(TMDs). While the properties of the 3D bulk materi¬ 
als have been well known for decades, the possibility of 
thinning them down towards the monolayer has given 
rise to an entirely new research area. Their structural 
formula MX 2 (with M being a transition metal and X a 
chalcogenide) comprises metals, semimetals, semiconduc¬ 
tors and superconductors [T] . What they have in common 
is that, with approaching the 2D limit, a whole world of 
intriguing properties such as extraordinarily large exci- 
ton binding energies [2] , robust valley polarization [3] and 
large spin orbit splitting [3] opens up. These properties 
pave the way for potential use of TMDs for applications 
in digital electronics and optoelectronics [3, in energy 
conversion and storage as well as in spintronics [7]. 
Most tungsten and molybdenum based TMDs exhibit a 
transition from an indirect to a direct bandgap semicon¬ 
ductor when being thinned down to the monolayer, re¬ 
sulting in high intensity photoluminescence [HHini- 
Raman spectroscopy is one of the most powerful tools in 
characterizing nanomaterials. For few-layer (FL-) TMDs 
such as FL-WS 2 it allows for instance to exactly deter¬ 
mine the number of layers. To be able to extract all the 
information that the measured Raman spectra have to of¬ 
fer, it is of primary importance to have a common ground 
on which to categorize the different Raman modes as a 
function of the number of layers. There is a large number 
of publications on monolayer TMDs mHH] and the re¬ 
spective bulk materials [HIISIIIS], but only few articles 
concentrate on the transition from monolayer to bulk, 
e.g. the evolvement of the Raman signatures with the 
number of layers. In these articles moreover, many use 
the symmetry of the bulk to assign the Raman features 
[iniinHil]; only recently there have been some reports 


that take into account the different symmetries of even 
number of layers (even N) and odd number of layers (odd 
N) TMDs [^nidl]- However they only focus on Raman 
modes that are allowed in first-order scattering. 

In this work we will show that (i) it is important to dis¬ 
tinguish between even and odd number of layers and (ii) 
that when the excitation energy is in resonance with the 
first optical transition of the investigated material it is 
necessary to consider the full set of phonons. We study 
the splitting of the monolayer out-of-plane mode in 
FL-WS 2 in particular and are able to observe layer de¬ 
pendent Raman signatures comprising Raman active and 
inactive modes. They allow for an easy identification of 
the number of layers via Raman spectroscopy. More¬ 
over a general systematic behavior for the splitting of 
the monolayer Raman modes is proven experimentally 
for the first time and supported by DFT calculations. 
The findings described in this work should be expand¬ 
able to other FL-TMDs. 


Experimental 

The samples are prepared from a bulk WS 2 crystal (hq 
graphene, Netherlands - Groningen) using the mechan¬ 
ical exfoliation technique. To enhance the optical con¬ 
trast, the crystals are exfoliated onto a 90 nm Si02/Si 
wafer. Raman measurements on WS 2 samples were done 
at room temperature in backscattering geometry with a 
Horiba Jobin Yvon LabRAM HR spectrometer using a 
confocal setup with a lOOx objective and excitation wave¬ 
lengths of 457 nm and 633 nm. To avoid sample heating 
the laser power was kept below a maximum of 120/iW; 
a 1800 lines mm“^ grid was used to ensure high spectral 
resolution of around lcm“^. First, Raman spectra were 
taken with the setup described above and calibrated with 
Neon lines. In a second step the measurements were re- 
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F-point were calculated in the frame of density functional 
(perturbation) theory on the level of the local density 
approximation (LDA) as implemented in the CASTEP 
code [53]. We treated the W(4d,5s) and the S(3s,3p) 
states as valence electrons using normconserving pseu¬ 
dopotentials with a cutoff energy of 800 eV. All recip¬ 
rocal space integrations were performed by a discrete 
fc-point sampling of 18x18x1 fc-points in the Brillouin 
zone. Starting from fully symmetric model geometries 
of one to five layers of AB-stacked WS 2 , we fully op¬ 
timized the lattice constants and atomic positions un¬ 
til the residual forces between atoms were smaller than 
0.01 eV/A and the stresses on the cell boundaries were 
smaller than 2.5xl0“^GPa. The obtained in-plane lat¬ 
tice constants slightly increased with layer number from 
a value of 3.141 A for IL-WS 2 to 3.143 A for 5 L-WS 2 , 
in good agreement with the experimental in-plane lattice 
constant of 3.15 A in bulk WS 2 |24|. Interactions of the 
sheet with residual periodic images due to the 3D bound¬ 
ary conditions were minimized by maintaining a vacuum 
layer of at least 20 A. 


FIG. 1: (a) Optical micrograph of a WS 2 flake with selected 
regions inspected by AFM [shown in (b) and (c)].(b) and (c) 
AFM images of the selected regions with indication of layer 
numbers and lines showing where AFM height profiles were 
taken from, (d) AFM height profiles going from the substrate 
onto the flakes for different positions on the flake. Step heights 
in units of the observed interplanar spacing of ~0.7nm are 
indicated. 

peated in subpixel (6sp) mode. There, each spectrum is 
taken a couple of times, while each time the spectrometer 
is shifted by a step size which is smaller than a pixel value. 
This technique does not increase the spectral resolution, 
however, by providing more data points per wavenum¬ 
ber, it reduces the signal to noise ratio significantly. The 
spectra aquired in subpixel mode were then shifted in fre¬ 
quency to match the calibrated Raman spectra obtained 
in the regular single-window mode. For better compari¬ 
son, in Figs. 2, 3, 4, the spectra were normalized to the 
intensity of the out-of-plane mode. 

Atomic force microscopy (AFM) images were acquired 
using a Park Systems XE-100 setup with commercial sil¬ 
icon tips in tapping mode configuration. Images were 
taken with 256x256 px resolution (4x4 fim). An exem¬ 
plary AFM analysis of a FL-WS 2 sample is shown in 
Fig. 1. Step heights between subsequent layer numbers 
were typically around 0.7 nm, close to the experimental 
value of the interplanar spacing of WS 2 layers |55|. Down 
to the monolayer, an offset of typically around 1 nm be¬ 
tween substrate and sample was observed, which is prob¬ 
ably due to the presence of adsorbates in between the 
substrate and the sample and was taken into account in 
the analysis presented in Fig. 1. 

The phonon frequencies of monolayer and FL-WS 2 at the 


Results 

WS 2 , like M 0 S 2 , crystallizes in the 2H trigonal pris¬ 
matic structure where the tungsten atoms are sand¬ 
wiched between two layers of sulfur atoms. Intralayer 
bonds are of covalent nature, the interlayer interaction is 
governed by weak van-der-Waals forces. Bulk WS 2 has 
the Dgft, point group; the 6 atoms per unit cell result in 18 
phonon modes at the T-point of the hexagonal Brillouin 
zone [55] : 

DqIi ■ r = Aig + 2A2„ -I- -I- 2i32g + Eig + 2Eiu -b E 2 U + 

2E2g 

In the monolayer and an odd number of layers (odd A), 
the symmetry is reduced to the point group. There¬ 
fore, odd N WS 2 do not have a center of inversion. The 
F-point phonon modes transform according to the follow¬ 
ing irreducible representation: 

D3h--'^= ^ {A[ + A" + E' + E") + A" + E', 

N = 1,3,5,... 

For even number of layers (even A), WS 2 possesses a cen¬ 
ter of inversion, the symmetry is described by the point 
group Da^: 

Dzd '■ F = ^(Alg + Eg + A 2 U + Eu), 

A = 2,4,6,... 

Let us now first consider the case where the excitation 
wavelength is far from resonance. Figure 2 (a) shows Ra¬ 
man spectra of one to five WS 2 layers and the bulk mate¬ 
rial taken with an excitation wavelength of 457 nm. The 
two main Raman modes are the A)^ and Aig mode around 
420 cm“^ for odd and even A, respectively, and the E' and 
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FIG. 2: (a) Raman spectra of FL WS 2 from one (IL) to five (5L) layers and the bulk material (NL). Excitation wavelength 
is 457 nm. Spectra are normalized to the intensity of the A mode at 420 cm“^. (b) Evolution of the frequency of the main E 
and A mode with the number of layers. The dashed line marks the increasing frequency difference between the modes with 
the number of layers, (c) Strong photoluminescence is observed for monolayer WS 2 at the direct excitonic transition energy 
around 625 nm. The inset shows a magnified image of the same spectrum in the region of the main IL-WS 2 Raman peaks. 


mode around 355 cm“^ for odd N and even N, respec¬ 
tively. A clear upshift with increasing number of layers 
is seen for the out-of-plane A\/Aig mode, whereas the 
in-plane E'/E^g mode slightly softens. This has been ob¬ 
served before for WS 2 [l9l|26] as well as for other TMDs 
[iniiniEilizH]. The stiffening of the out-of-plane A mode 
is explained by the increasing interlayer interaction and 
the subsequent rise in restoring forces on the atoms with 
the number of layers m- The same should hold for the 
E mode, although to a lesser extent, as the atoms move 
in-plane and the influence of interlayer interaction is thus 
expected to be smaller. However, the opposite trend is 
observed and has been attributed to dielectric screen¬ 
ing of long range Coulomb interactions [T^. Figure 2 
(b) depicts the change in frequency of the out-of-plane 
A and in-plane E mode with the number of layers and 
the frequency difference between the two modes, which 
increases from 60.4cm“^ to 65cm“^ from the monolayer 
to the bulk. 

If the exciting light is far from resonance, the spectra 
are dominated by Raman modes allowed in first-order 
scattering, as has been reported previously for WS 2 nan¬ 
otubes [55]. If we focus on the out-of-plane A modes, it is 
evident that the mono-, bilayer and bulk spectra show a 
single peak, whereas for three and more layers there is at 
least another mode appearing as a low-energy shoulder 
of the dominant Raman feature. This apparent splitting 
of the out-of-plane mode is due to the fact that for mono- 
layer and bulk WS 2 there is only one Raman active 
and Aig mode, respectively; for few layers starting with 
three layers, more than one Raman mode becomes al¬ 
lowed, see also Ref. [2T] . 

To further explore these new Raman modes, which are 
reported here for the first time in FL-WS 2 , we analyze 
the Raman spectra of the same samples taken under the 


resonance condition. Figure 2 (c) shows a photolumi¬ 
nescence spectrum of monolayer WS 2 at 457 nm excita¬ 
tion wavelength. The photoluminescence signal is more 
than two orders of magnitude larger than the Raman 
signal and has maximum intensity around 625 nm. From 
previous experiments it is known that the first optical 
transition energy of bulk WS 2 is constituted by the A 
exciton around 633 nm |29j . Therefore, with 633 nm ex¬ 
citation wavelength, we are close to the A excitonic res¬ 
onance for mono- and few layer WS 2 . In the remainder 
of this work we will focus on the out-of-plane A mode 
around 420 cm“^. Of the four Raman modes allowed 
in bulk WS 2 Deh symmetry, the E^g is too low in fre¬ 
quency to be observed here, the Fig mode is not allowed 
in backscattering geometry, and the E^g mode around 
350 cm“^ overlaps with a second order mode, which dom¬ 
inates the spectra in resonance m- As we will show be¬ 
low, the out-of-plane mode is clearly separated from other 
Raman features and shows sidebands that can readily be 
explained by the respective symmetries of even or odd 
number of layers. 

Figure 3 shows the region of the out-of-plane Raman 
mode of few-layer WS 2 taken with 633 nm excitation 
wavelength. Similarly to the spectra shown at 457 nm 
excitation wavelength, the upshift of the Raman mode 
with the number of layers is evident. More importantly 
though, there are striking differences in the shape of the 
Raman mode. For FL-WS 2 the structure of the out-of- 
plane mode is complex. In contrast to the spectra de¬ 
scribed above, already the bilayer spectrum shows more 
than one component. More and more sidebands arise for 
increasing number of layers but they get weaker for n > 5 
and vanish for very thick flakes (n = TV). This again un¬ 
derlines the special role played by few-layer samples: for 
n = 1 and n = N the symmetry of WS 2 allows only one 
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FIG. 3: Experimental resonance Raman spectra (excitation 
wavelength 633 nm) of the A mode region in FL-WS 2 from 
the bilayer (2L) to five layers (5L) and the bulk (NL). The 
monolayer spectrum (IL) is only shown with an excitation 
wavelength of 457 nm, as the Raman features are dominated 
by strong photoluminescence at 633 nm excitation spectrum 
(not shown). For FL-WS 2 , spectra of the same region taken 
with 457 nm excitation wavelength are overlayed and depicted 
with dashed lines. Spectra are normalized to the main Raman 
peak and are offset for clarity. 

Aj (n = 1) and Aig (bulk) Raman mode. The compari¬ 
son with the out-of-resonance spectra (457 nm excitation 
wavelength) illustrates that, even though shoulders of the 
main Raman peak are observed at 457 nm excitation as 
well, an increased number of well pronounced sidebands 
to the main Raman peak appears mainly for Raman mea¬ 
surements in resonance with the A exciton. We fit the 
spectra with Lorentzian profiles, see Fig. 4; for clarity 
spectra of even and odd N are shown in separate graphs. 
In Fig. 4 (a) the bilayer Aig mode spectrum possesses 
a low energy shoulder that has not previously been ob¬ 
served. Its appearance is surprising since the only ex¬ 
pected Raman active vibration is the Aig mode, where 
both layers vibrate in-phase according to the monolayer 
A'l mode. We will show below that this second peak is in¬ 
deed the infrared active A 2 U mode, where the two layers 
vibrate out-of-phase. In the four-layer (4L) spectrum the 
dominant Aig mode shifts up with respect to the bilayer, 
and a prominent shoulder is seen at almost the same fre¬ 


quency as the shoulder in the bilayer spectrum. The 4L 
spectrum is best fitted with four Lorentzians to account 
for the Raman intensity between the two stronger Ra¬ 
man features. In the spectra of odd N WS 2 , a similar 
pattern is revealed. In Fig. 4 (b) the trilayer (3L) and 
five-layer (5L) spectra are shown. The former consists 
of the expected A'l peak as the most significant contri¬ 
bution and a pronounced low-energy shoulder. A third 
Lorentzian fits the plateau between the main peaks. The 
3L spectrum cannot be properly fitted with only two 
Lorentzians. From the experience gained from the spec¬ 
tra investigated above, the five-layer spectrum is fitted 
with five Lorentzians, two of which fill up the region be¬ 
tween the main A'l peak and the two low-energy shoul¬ 
ders. From the spectra for layer numbers of n = I to 
n = 5 it thus seems that there are always N compo¬ 
nents to the out-of-plane A mode, where A is the num¬ 
ber of layers. We have observed sidebands on the lower 
energy side of the dominant Raman peak also for higher 
layer numbers as exemplary shown for the spectrum of an 
around 10 layer thick flake. But they are rather weak and 
cannot be fitted following the pattern described above. 

DISCUSSION 

Recently, careful analysis of few-layer TMDs has led 
to the observation of Raman modes that are neither 
seen in the bulk nor in the monolayer m- Some of 
them appear as shoulders to Raman modes that are 
allowed in first order for bulk and monolayer, like the 
bulk Aig mode discussed here. Others are Raman 
inactive or not allowed in backscattering geometry in 
bulk and monolayer, like the bulk B 2 g and Eig modes 
[B HD 1301 HI]. For the first case, to the best of our 
knowledge, there is only one article that explicitly 
shows a splitting of the first-order Aig mode with the 
number of layers in few-layer MoSe 2 [10] • The 3L and 
4L samples show two components, for the 5L sample a 
third component is seen. For WSe 2 the overlap of the 
bulk E 2 g and Aig mode makes an observation of such 
shoulders impossible jOD ED] , and for the case of MoTe 2 
the intensity of the Aig mode appears to be too small 
to resolve a multipeak structure [05] . For M 0 S 2 there is 
only little literature on resonance Raman spectra, and 
despite some asymmetry in the shape of the out-of-plane 
Raman mode, a splitting similar to the one investigated 
in this work, is not observed [DDES]. Terrones et al. [21] 
calculate the optical phonons for a number of FL-TMDs, 
among them M 0 S 2 and WS 2 , but do not discuss phonons 
other than the Raman active ones. 

In order to have a theoretical background to the ex¬ 
perimentally observed appearance of more than just 
the Raman active vibrational modes in the spectra of 
FL-WS 2 , we have done calculations employing density 
functional theory (DFT). A better insight into the 
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FIG. 4: Experimental resonance Raman spectra of FL-WS 2 . The Lorentzian fit curves are shown as well as the symmetry 
attributed to the individual modes. In general, there are N components for a N layer spectrum, (a) For even number of layers, 
Raman active Aig modes alternate with infrared active A 2 U modes, (b) For odd number of layers, Raman active A'l modes 
alternate with infrared A 2 modes. For both, even and odd N, the Raman active modes are more intense than the rather weak 
infrared active modes. 


atomic displacements corresponding to the phonon 
modes of FL-TMDs is given in Fig. 5, where schematic 
drawings of all possible vibrations evolving from the 
monolayer mode in WS 2 from n= 1 — 5 are depicted, 
based on the DFT calculations. While there is only one 
possibility in monolayer WS 2 for the sulfur atoms to 
vibrate against each other with a fixed tungsten atom 
in between, a splitting of this mode occurs for bilayer 
WS 2 . Since bilayer WS 2 possesses a center of inversion, 
there is a Raman active Aig mode, where the two layers 
vibrate in phase, and an infrared active A 2 „ mode, 
where the two layers vibrate out of phase. As the latter 
is not Raman active it is not seen in Raman spectra 
taken with excitation wavelengths far from resonance 
[Fig. 2 (a)]. However, it is observed for the resonance 
Raman spectrum (Figs. 3 and 4), albeit with weaker 
intensity than the dominant Aig mode. Several possible 
reasons for this unusual behavior are discussed below. 
As the two layers interact more strongly for the in-phase 
vibration, the Aig mode has a slightly higher frequency 
than its infrared active counterpart. For 4 L-WS 2 , 
each of the two bilayer modes again splits up into a 
Raman-active Aig and a infrared active A 2 U mode. The 
spectrum is still governed by the in-phase vibration of 
all four layers (Aig) but there is a second Raman active 
Aig mode that has the outer layers vibrating out of 
phase with the inner ones, thus retaining the inversion 
symmetry of the overall structure. It is interesting to 
note that this lower lying Aig mode in the 4L spectrum 
has almost the same frequency as the infrared active 
mode in the bilayer, a pattern that will also be observed 


for odd N, see below. In addition, we identify the two 
small shoulders on the lower-frequency side of the two 
Raman active modes with the A 2 „ modes [Fig. 4 (a)], 
the lowest lying with neighboring layers vibrating out of 
phase, the other one with the two upper layers vibrating 
out of phase with the two lower layers. 

In odd N WS 2 there is obviously again the possibility 
of all layers vibrating in-phase and out-of-phase. In 
contrast to even N WS 2 , where the pure out of phase 
vibration is not Raman active, for odd N both, in- 
and out-of-phase vibrations, are Raman active and 
possess A'l symmetry. For trilayer WS 2 , the atomic 
displacement vectors of these two modes are shown on 
the right and left side of the third panel of Fig. 5. In 
the spectrum depicted in Fig. 4 (b), the lower lying A'l 
mode accounts for the strong shoulder at approximately 
416cm“^ of the main A'l mode (in-phase vibration). 
In between the modes a plateau is evident that is not 
accounted for if the spectrum is only fitted with two 
Lorentzians. The origin of the plateau is attributed to an 
infrared active A 2 mode, with the middle layer fixed and 
the sulfur atoms of the top and bottom layer vibrating 
out-of-phase (see Fig. 5, third panel, middle). The same 
approach can now be used for the analysis of the five 
layer WS 2 spectrum. Two shoulders to the main A'^^ 
peak can be identified and, following the pattern of Fig. 
5, attributed to another two Raman active A'l modes. 
In between the Raman active modes two very weak 
features belong to infrared active A 2 vibrations. Coming 
from the trilayer WS 2 , the two modes with the highest 
frequency in the five-layer material can be imagined as 
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FIG. 5: Schematic drawing for all possible vibrational modes 
in the out-of-plane mode region of FL-WS 2 (Raman and in¬ 
frared) for one (IL) to five (5L) layers together with the sym¬ 
metry assignments taken from our DFT calculations. The 
displacement patterns are ordered with increasing frequency 
from left to right. Taking into account the full set of possi¬ 
ble vibrational patterns helps to study the splitting of vibra¬ 
tional modes with increasing layer number and to attribute 
the modes to the features seen in the Raman spectra. 


stemming from a splitting of the main Raman active 
A'l mode in trilayer. The same can be said about the 
lower frequency Raman active mode in the trilayer that 
splits up into the low frequency A'l and A 2 mode in five 
layers. The only infrared active A 2 mode in the trilayer 
spectrum changes symmetry in the hve layer and has 
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TABLE I: Experimental Raman frequencies (in cm“^) of all 
vibrational modes observed in FL-WS 2 excited with 633 nm. 
For N layers the A'l mode of the monolayer splits up into 
N components. The point group of FL-WS 2 is Dsh for odd 
number of layers and D^d for even number of layers. The vi¬ 
brational modes for every number of layers alternate between 
Raman and infrared active. 


the two outermost layers vibrating out-of-phase with the 
three inner ones. 

To summarize, for analyzing the Raman spectra of one 
to five layers of WS 2 , it is necessary to take into account 
the full set of phonon modes. For the evolvement of the 
monolayer A'j^ mode analysed in detail here, the number 
of out-of-plane modes matches the number of layers. 
Raman active modes alternate with modes that are 
infrared active and are seen in the spectra of FL-WS 2 for 
the first time. The pattern observed for the out-of-plane 
mode should also be valid for all other vibrational modes 
in FL-WS 2 and other FL-TMD materials. For even N, 
the monolayer mode splits up into ^ Raman active and 
^ infrared active modes (the in-plane modes are always 
doubly degenerate). For odd N, the monolayer mode 
evolves into Raman active and infrared active 
components. In total, in N layers, each of the monolayer 
phonons splits up into N phonon modes EHIEI]. 

Table 1 lists the experimentally obtained Raman fre¬ 
quencies of all out-of-plane vibrations derived from the 
Aj mode at 418.8 cm“^ for few-layer WS 2 . Where two 
or more samples with the same number of layers were 
measured, the average value is given in Table 1. In all 
these cases deviations from the given frequencies are less 
than 0.3 cm“^. Again, the table illustrates that Raman 
and infrared active modes are alternating irrespective of 
the layer number. Additionally, not only the main Ra¬ 
man active component with A'j^/Aig symmetry for odd 
and even N appears to exhibit increased frequency with 
increasing number of layers but also the other compo¬ 
nents (2nd row and below in Table 1) follow the same 
trend. This is illustrated in Fig. 6 (a), where the tabu¬ 
lated frequencies are plotted against the number of lay¬ 
ers. The main Raman peak in the spectra shown in Figs. 
2-4 is seen stiffening in frequency from the monolayer 
to five layers (red circles, connected by a dashed line to 
guide the eye). Starting with the bilayer an infrared ac¬ 
tive mode comes into play that also appears for higher 
layer numbers. It also increases in frequency, thus fol¬ 
lowing the behavior of the main Raman component due 
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FIG. 6: a) Raman frequencies of the A type modes in FL-WS 2 . Red circles mark the Raman active modes, blue squares denote 
the infrared active modes also seen in the Raman spectra. Dashed lines are used to guide the eye. Starting from one layer the 
splitting of the modes for FL-WS 2 produces a fan-like shape of the possible vibrational frequencies, b) Theoretical calculations 
with DFT reproduce the observed trends. Raman active modes are denoted with open red circles, infrared active modes are 
marked with open blue squares. 


to increased force constants with increasing number of 
layers (blue squares, connected by a dashed line to guide 
the eye). The same is seen for a second Raman active 
feature starting with three layers and another infrared 
active mode starting with four layers. Interestingly, the 
frequencies observed in the mono- , bi- and trilayer are al¬ 
most exactly repeated when the layer number is increased 
by two, underlining the close relation of the out-of-plane 
modes even though the symmetry and Raman/infrared 
activity changes with the layer number. In contrast, the 
position of the lowest frequency mode stays almost con¬ 
stant from three layer onwards. This mode always has 
neighboring layers moving out-of-phase, but neighboring 
sulfur atoms from adjacent layers moving in-phase. As a 
result, the nearest neighbor force constants determining 
the frequency of this mode will not change significantly 
for larger number of layers. In Fig. 6 (b), the frequen¬ 
cies calculated with DFT are plotted against the layer 
number. Despite a slight overestimation of absolute fre¬ 
quencies and a smaller magnitude of the splitting of the 
modes, the experimental results are well reproduced. As 
a whole, the splitting of the out-of-plane monolayer A)^ 
mode in FL-WS 2 results in a fan-like shape showing simi¬ 
larity to Fig. 5 in Ref. |2D]. Zhang et aL ^0] investigated 
the evolvement of the low-frequency rigid layer C (dis¬ 
placement along the c-direction) and LB (layer breath¬ 
ing) modes in FL-M 0 S 2 , a first report on these interlayer 


modes in few-layer M 0 S 2 can be found in Ref. [34]). With 
the support of a simple atomic chain model the authors 
found a behavior similar to the one described here for 
the out-of-plane vibration in WS 2 . Obviously, rigid layer 
vibrations only appear starting from the bilayer. Strictly 
speaking, the difference to all other optical modes in FL- 
TMDs is that there are consequently only A-1 possible 
vibrations (where N denotes the layer number) for the 
rigid layer modes (E^g and B 2 g symmetry in bulk M 0 S 2 ). 
Of course, the “N modes for N layer rule” is restored 
again if one adds the acoustic modes. The acoustic E' 
and A 2 modes of the monolayer are the origin of the low 
frequency vibrations in FL-TMDs and will split up into 
N components for N layers (among them an A and an E 
type mode with zero frequency). 

The behavior described above should in principle be ob¬ 
servable in other 2H-TMDs as well. Terrones et al. [21] 
predict an increased splitting of the Raman active compo¬ 
nents in the order WSe 2 , MoSe 2 , WS 2 and M 0 S 2 . What 
distinguishes WS 2 from other prominent 2H-TMDs like 
WSe 2 and M 0 S 2 is that the bulk Aig and E^g modes are 
well separated in energy and that no second-order Ra¬ 
man features overlap with the Aig mode, thus making 
it easier to resolve the splitting of the out-of-plane mode 
into Raman and infrared active components. More im¬ 
portantly, measuring in resonance with the optical tran¬ 
sition appears to be a necessary condition to observe the 





full set of vibrational modes. There is a lack of studies 
on Raman spectra of other FL-2H-TMDs measured under 
resonance condition; often the corresponding excitation 
wavelengths are avoided because in the monolayer case 
the Raman features are obscured by strong photolumi¬ 
nescence signal. In FL-WS 2 in particular, the character¬ 
istic shape of the out-of-plane mode in resonance Raman 
spectra can be used as a fingerprint region to unambigu¬ 
ously identify the number of layers. 

For WS 2 measured under the resonance excitation, even 
in the bulk material, the Aig out-of-plane mode is ac¬ 
companied by a small shoulder that is attributed to the 
silent Bln mode [121 [161 [29] . It had been previously at¬ 
tributed to a LA(K) + TA(K) combination mode |35|, but 
the participation of two phonon modes involving acous¬ 
tic phonons in this frequency region seems unlikely, espe¬ 
cially in light of a more recent DFT calculations showing 
that nowhere in the Brillouin zone the acousting phonons 
reach values above 200 cm“^ [12]. The in-phase and out- 
of-phase vibrational pattern of the Aig/Bi„ pair in bulk 
WS 2 finds its counterpart in the variety of in- and out- 
of-phase vibrations observed in FL-WS 2 . In WS 2 nano¬ 
materials, the Bitj mode gains in Raman intensity and 
its evolvement can be followed in WS 2 nanomaterials 
under pressure |29L I36j , in different layer orientation in 
thin films m and in different diameter WS 2 nanotubes 
[221 [3H]- In these cases, a strong resonance behavior 
is observed as well, much like in FL-WS 2 : the modes 
not allowed in a first-order Raman process appear most 
strongly when the excitation energy is in or close to res¬ 
onance with the optical transitions. In an earlier work 
on WS 2 nanotubes [22], we found that the curvature- 
induced strain and the resultant crystal symmetry dis¬ 
tortion was to be held accountable for the activation of 
silent modes. Here, in quasi 2D materials, the situation 
in absence of curvature effects is different; strain due to 
substrate-sample interaction is assumed to only play a 
very minor role in Raman spectroscopy on supported FL- 
TMDs [22]. 

Instead, a closer look at the nature of the excitonic transi¬ 
tion leading to the resonantly enhanced Raman intensity 
in the spectra of FL-WS 2 can provide means to elucidate 
the appearance of the infrared-active components of the 
Aj/Aig Raman mode around 420 cm“^. For 633 nm ex¬ 
citation wavelength, the phonons are coupling to the A 
exciton situated at the K point of the Brillouin zone. 
Owing to the structural relationship between strong in¬ 
tralayer bonding in two dimensions and weak interlayer 
interaction in the third dimension, like many properties 
of layered TMD systems, the wavefunctions of the exci- 
tons are expected to be very anisotropic. A recent work 
on the orientation of luminescent excitons in FL-M 0 S 2 , 
isostructural to FL-WS 2 , reveals them being confined en¬ 
tirely in-plane without significant expansion in the stack¬ 
ing direction of individual layers |39] . This was further 
substantiated by density-functional theory (DFT) calcu¬ 


lations in the local-density approximation (LDA) that 
explicitely showed the wavefunction of the A exciton in 
multilayer M 0 S 2 to be spread out over a large area in 
two dimensions but with negligible density in neighbour¬ 
ing layers [JD] • For the resonant Raman process discussed 
here this means that even for a layer number of more than 
one, the phonon couples to an A exciton localized primar¬ 
ily in one of the layers. If the N layers in FL-WS 2 are 
to be treated approximately as N individual monolayers 
for the specific case of the A excitonic resonant Raman 
process, there are N Raman allowed Raman modes to be 
expected in the region of the monolayer A)^ Raman mode 
with similar Raman intensities. They are still split in 
frequency due to interlayer interaction. Here, all N com¬ 
ponents are identified in the Raman measurements pre¬ 
sented in this article, but the infrared-active components 
- speaking from the TV-layer symmetry point of view - are 
always weaker in intensity than the Raman active com¬ 
ponents. Thus we conclude that Raman selection rules 
are at least weakened but not completely broken. This is 
supported by the fact that the Raman active components 
are gaining in intensity relative to the main Raman peak 
but still appear as shoulders rather than as individual 
peaks. 

Far from the A excitonic resonance, the few-layer WS 2 
can not be treated as N individual monolayers and the 
Raman selection rules following from the few-layer sym¬ 
metry strictly apply. This is also in agreement with re¬ 
cent findings on newly observed Raman modes in FL- 
M 0 S 2 in resonance with the C exciton [31] . 


Conclusion 

In summary, we have shown experimentally as well as 
theoretically that the out-of-plane A^ mode of the WS 2 
monolayer splits up in the few layer regime into N com¬ 
ponents for N layers. Despite the fact that only A of 
them for even number of layers and for odd num¬ 
ber of layers are Raman active, the full set of phonon 
modes is observed when the laser excitation energy is 
close to the A excitonic transition energy. A possible ex¬ 
planation for this unusual behavior is presented by tak¬ 
ing into account the in-plane orientation of the A exciton 
wavefunction involved in the resonant Raman scattering 
process. The stiffening of the main out-of-plane phonon 
mode with N is also followed by all other components 
successively added with increasing number of layers. By 
resonant Raman scattering measurements one can con¬ 
clusively identify the number of layers in a specific sample 
by simply counting the number of components of the out- 
of-plane A mode. The detailed analysis of the evolvement 
of the Aj mode of monolayer WS 2 presented here should 
in principle be applicable to (i) all other Raman modes 
of (ii) all layered materials in the few-layer regime. 


9 


Acknowledgements 

This work was supported by the Deutsche Forschungs- 
gemeinschaft (DFG) in the Priority Programme SPP 
1459 “Graphene” and by the European Research Council 
(ERG) through grant number 259286. We thank Lukas 
Madaufi for his support with sample preparation. 


* Corresponding author, nistaiger@physik.tu-berlin.de 

[1] M. Chhowalla et al., Nat. Chem. 5, 263 (2013). 

[2] Z. Ye et al., Nature 513, 214 (2014). 

[3] B. Zhu, H. Zeng, J. Dai, Z. Gong, and X. Cui, Proc. 
Natl. Acad. Sci. Ill, 11606 (2014). 

[4] A. Ramasubramaniam, Phys. Rev. B 86, 115409 (2012). 

[5] Q. Wang, K. Kalantar-Zadeh, A. Kis, J. Coleman, and 
M. S. Strano, Nat. Nanotech. 7, 699 (2012). 

[6] H. Wang, H. Feng, and J. Li, Small 10, 2165 (2014). 

[7] B. Zhu, H. Zeng, J. Dai, and X. Cui, Adv. Mat. 26, 5504 
(2014). 

[8] A. Splendiani et al., Nano Lett. 10, 1271 (2010). 

[9] K. F. Mak and C. Lee and J. Hone and J. Shan and 
T. F. Heinz, Phys. Rev. Lett. 105, 136805 (2010). 

[10] P. Tonndorf et al., Optics. Expr. 4, 4908 (2013). 

[11] H. Matte et al., Angew. Chem. Int. Ed. 122, 4153 (2010). 

[12] A. Molina-Sanchez and L. Wirtz, Phys. Rev. B 84, 
155413 (2011). 

[13] M. Thripuranthaka and J. Dattatray, ACS Appl. Mater. 
Interfaces 6, 1158 (2014). 

[14] A. A. Mitioglu et ah, Phys. Rev. B 89, 245442 (2014). 

[15] T. Wieting and J. Verble, Phys. Rev. B 3, 4286 (1971). 

[16] C. Sourisseau, M. Fouassier, M. Alba, A. Ghorayeb, and 


O. Gorochov, Mater. Sci. Eng., B 3, 119 (1989). 

[17] C. Lee et al., AGS Nano 4, 2695 (2010). 

[18] G. Plechingeret al., Appl. Phys. Lett 101, 101906 (2012). 

[19] A. Berkdemir et al.. Sc. Rep. 3, 1755 (2013). 

[20] X. Zhang et al., Phys. Rev. B 87, 115413 (2013). 

[21] H. Terrones et ah. Sc. Rep. 4, 4215 (2014). 

[22] J. Wilson and A. Yoffe, Adv. Phys. 18, 193 (1969). 

[23] S. Glark et al., Z. Kristallogr. 220, 567 (2005). 

[24] W. Schutte, J. D. Boer, and F. Jellinek, J. Solid. State 
Chem. 70, 207 (1987). 

[25] J. Verble and T. Wieting, Phys. Rev. Lett. 25, 362 
(1970). 

[26] W. Zhao et al.. Nanoscale 5, 9677 (2013). 

[27] D. Late, B. Liu, H. Matte, C. Rao, and V. Dravid, Adv. 
Funct. Mater. 22, 1894 (2012). 

[28] M. Yamamoto et al., ACS Nano 8, 3895 (2014). 

[29] M. Staiger et al., Phys. Rev. B 86, 165423 (2012). 

[30] X. Luo, Y. Zhao, J. Zhang, Q. Xiong, and S. Y Quek, 
Phys. Rev. B 88, 075320 (2013). 

[31] N. Scheuschner, R. Gillen, M. Staiger, and J. Maultzsch, 
in preparation (2015). 

[32] N. Scheuschner, O. Ochedowski, M. Schleberger, and 
J. Maultzsch, Phys. Status Solid! B 249, 2644 (2012). 

[33] H. Li et al., Adv. Funct. Mater. 22, 1385 (2012). 

[34] H. Zeng et al., Phys. Rev. B 86, 241301(R) (2012). 

[35] C. Sourisseau et al., Ghem. Phys. 150, 281 (1991). 

[36] S. Yu et al., J. Phys.: Gondens. Matter 19, 425228 
(2007). 

[37] J. Ghung, Z. Dai, K. Adib, and F. Ohuchi, Thin Solid 
Films 335, 106 (1998). 

[38] M. Krause et al., ChemPhysChem 10, 2221 (2009). 

[39] J. A. Schuller et al., Nat. Nanotechnol. 8, 271 (2013). 

[40] A. Molina-Sanchez et al., Phys. Rev. B 88, 045412 
(2013). 



